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THE FOURTH-ORDER TYPE LINEAR 
ORDINARY DIFFERENTIAL EQUATIONS 

W.N. EVERITT, D. J. SMITH, AND M. VAN HOEIJ 

Abstract. This note reports on the recent advancements in the search for exphcit repre- 
sentation, in classical special functions, of the solutions of the fourth-order linear ordinary 
differential equations named Bessel-type, Jacobi-type, Laguerre-type, Legendre-type. 
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1. Introduction 

The four fourth-order type linear ordinary differential equations are named after the four 
classical second-order equations of Bessel, Jacobi, Laguerre, Legendre. 

The three fourth-order equations Jacobi-type, Laguerre-type, Legendre-type, and their 
associated orthogonal polynomials, were first defined by H.L. Krall in 1940, see ^U] and |llj . 
and later studied in detail by A.M. Krall in 1981, see jH], and by Koornwinder in 1984, see 



The structured definition of the general-even order Bessel-type special functions is depen- 
dent upon the Jacobi and Laguerre classical orthogonal polynomials, and the Jacobi-type 
and Laguerre-type orthogonal polynomials; see ,8J. 

The properties of the fourth-order Bessel-type functions have been studied in the papers 
0, and 0. 

This short paper reports on the recent progress that has been made in the search for 
representation of solutions of these four type differential equations in terms of classical special 
functions. 
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For the Bessel-type equation two independent solution were obtained, see and [3j, and 
then a complete set of four independent solutions in the papers [Hj and these solutions 
are dependent upon the classical Bessel functions Jr,Yr,Ir, for r = 0, 1. 

Subsequently, using the techniques developed in jTSI, solutions of the three Jacobi-type, 
Laguerre-type, Legendre-type equations have been obtained in terms of the solutions of the 
corresponding classical differential equations. 

In all these four cases the verification of these representations as solutions of the type 
differential equations, is a formal process using the computational methods in the computer 
program Maple. Copies of the corresponding Maple .mws files, for the Jacobi-type, Laguerre- 
type, Legendre-type equations, are available; see Section below. 

There are four subsequent sections devoted to reporting on the form of these solutions, 
for each of the type fourth-order differential equations. 

The four type differential equations are all written in the Lagrange symmetric (formally 
self-adjoint) form with spectral parameter A or A; if further properties of the solutions of 
these equations, involving the theory of special functions, are to be studied, then these 
spectral notations may have to be changed. 

2. BeSSEL-TYPE DIFFERENTIAL EQUATION 

The fourth-order Bessel-type differential equation takes the form 

(2.1) {xy'\x))" - ((9s"^ + m-^x)y\x))' = Axy{x) for all x e (0, oo) 

where M G (0, oo) is a positive parameter and A G C, the complex field, is a spectral 
parameter. The differential equation (j2.1|) is derived in the paper Section 1, (1.10a)]. The 
many spectral properties of this equation in the weighted Hilbert function space L'^{{0, oo);x) 
are studied in the papers |3j, |2j and [H]; for a collected account see the survey paper [31. 

Our knowledge of the special function solutions of the Bessel-type differential equation 
(j2.1|) is now more complete than at the time the paper [H] was written. However, the results 
in |Hl Section 1, (1.8a)], with a = 0, show that the function defined by 

(2.2) Jl'^ix) := [1 + M(A/2) Vo(Aa;) - 2M{X/2)\Xx)-^Ji{Xx) for all x G (0, oo), 

is a solution of the differential equation ()2.1|) . for all A G C, and hence for all A G C, and all 
M > 0. Here: 

{{) the parameter M > 
(m) the parameter A G C 

(iii) the spectral parameter A and the parameter M, in the equation (j2.1|) . and the pa- 
rameters M and A, in the definition ()2.2|) . are connected by the relationship 

(2.3) A = A(A, M) = A2(A2 + SM"^) for all A G C and all M > 

(iv) Jo and Ji are the classical Bessel functions (of the first kind), see [13 Chapter HI]. 
Similar arguments to the methods given in [S] show that the function defined by 

(2.4) Y!^'^\x) := [1 + M{X/2f]Yo{Xx) - 2M{X/2)\Xx)-^Yi{Xx) for all x G (0, oo). 
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is also an independent solution of the differential equation ()2.1|) . for all A G C, and hence for 
all A e C and all M > 0; here, again, Yq and Yi are classical Bessel functions (of the second 
kind), see jTHl Chapter III]. 

These earlier studies of the fourth-order differential equation (j2.ip failed to find any ex- 
phcit form of two linearly independent solutions, additional to the solutions J°' and Y^' . 
However, results of van Hoeij, see [13] and [Hj, using the computer algebra program Maple 
have yielded the required two additional solutions, here given the notations of /°'^ and 
K^^\ with explicit representation in terms of the classical modified Bessel functions Io,Kq 
and Ii,Ki. These two additional solutions are defined as follows, where as far as possible we 
have followed the notation used for the solutions J^'^''^ and Y^'^ , 

(i) given A G C, with arg(A) G [0, 27r), M G (0, oo) and using the principal value of ^/^, 
define 

(2.5) c = c(A, M) := \/A2 + and d = d{X, M) := 1 + M(A/2)2 

(a) define the solution, for all x G (0, oo), 
(2.6) 

I^'^ix) := -dioicx) + IcMx'^hicx) 

(2.7) := -[1 + M(A/2)Vo (xVX^ + SAf-i) + ^/M {2 + M{X/2y)x-^h (xVX^ + 8M-i) 

(iii) define the solution, for all x G (0, oo), 
(2.8) 

kI'^{x) := dKo{cx) + \cMx'^Ki{cx) 
(2.9) 

:= [1 + M(A/2)2]Ko (xVA2 + 8M-i) + ^ M [2 + M{X/2Y)x-^ Ki [x^X^ + 8M-i) . 

Remark 2.1. We have 

(1) The four linearly independent solutions Y^'^'^ , provide a basis for all 
solutions of the original differential equation 1)2.11) . subject to the (A, A) connection 
given in ()2.Hj) . 

(2) These four solutions are real-valued on their domain (0, oo) for all A G M. 

(3) The domain (0, oo) of the solutions and /°'*^ can be extended to the closed 
half-line [0, oo) with the properties 

jO'A^(O) = /°'^^(0) = 1 for all A G M and all M G (0, oo). 
3. LAGUERRE-TYPE differential EQUATION 

The Laguerre-type differential equation was discovered by H.L. Krall, see [Hlj and [TT], 
and subsequently studied by other authors, see [HI, [12] and [7|. 

The differential equation may be written in two forms; here the parameter A G (0, oo) and 
the spectral parameter A G C: 
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(1) The Frobenius form: 

, . W (a;) + (_2a;2 + 4a;)y{3) (a;) + (a;2 _ (^2A + 6)x)y"{x) 

^ ' +{{2A + 2)x - 2A)y'{x) = Xy{x) for all x e (0, oo). 

(2) The Lagrange symmetric form (formally self-adjoint form): 

.3 2) ii^^ exp{-x)y"{x))" - (((2A + 2)x + 2) exp(-x)y'(a;))' 

= Xexp{—x)y{x) for all x G (0, 00). 

The spectral properties of the equation (|3.2p are considered in the weighted Hilbert func- 
tion space L'^{{0, 00); exp(-)). 

The van Hoeij method of searching for factors and solutions of this differential equation 
depends on the use of the computer algebra program Maple; see ^| and 

Following the solutions of the fourth-order Bessel-type obtained in |HJ Section 2] , and then 
in |3] on using the methods of J3], it has proved possible to express all solutions of the 
equations ()3.1|) and ()3.2|) in terms of the confluent hypergeometric functions. 

The special functions involved are the Whittaker functions M^^z) and PVk,^(z) as defined 
in the compendium p!| Chapter 13, Sections 13.1.31 to 13.1.34]. 

Define the functions $ and ^ , for the range of variables given in jT| Chapter 13]: 

(3.3) $(K,/i;z) :=M,,^(z) 

(3.4) ^{k,^i-z):=W,,^{z). 
Given the parameters A e (0, 00) and A G C define 

(3.5) r(A, A) ■= V4A2 + 4A + 1 + 4A. 

We designate the four linearly independent solutions of ()3.1|) . equivalently ()3.2|) . again the 
parameters A G (0, 00) and A G C, by 

Lj.(A, A; x) for all x G (0, 00) and r = 1, 2, 3, 4. 

These solutions are defined as follows, where d/dx denotes differentiation with respect to 
the independent variable x, 

Li(A, A- x) := (i + ir(A, A)) exp(ix)<l'(-A - ir(A, A), 0; x) 
^^■^^ exp(ix)<l>(-A - |r(A, A), 0; x)) . 
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L2(A, A; x) := (| - |r(A, A)) x'^^^ exp(|a:)$(-A + |r(A, A), 0; x) 
^^'^^ {x-'/^ expilxM-A + ir(A, A), 0; x)) . 



L 
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L3(A, A; x) := (I + ir(A, A)) exp(ix)^(-A - |r(A, A), 0; x) 
^^'^^ exp(ix)^(-A - ir(A, A), 0; x)) . 
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(3.9) 



L4(A, A; x) 



(i - ir(A, A)) exp(ix)^(-A + ir(A, A), 0; x) 
exp(ix)^(-A + ir(A, A), 0; x)) . 



4. LEGENDRE-TYPE differential EQUATION 

The Legendre-type differential equation was discovered by H.L. Krall, see (TU] and [TT| . 
and subsequently studied by other authors, see [H], [12] and [7j. 

The differential equation may be written in two forms; here the parameter A G (0, oo) and 
the spectral parameter A G C: 

(1) The Frobenius form: 

^ ' ^ +8Axy'{x) = Xy{x) for all x G (-1, +1). 

(2) The Lagrange symmetric form (formally self-adjoint form): 

((1 - x'Yy"ix))" - ((8 + 4A(1 - x'))y'ix))' 
^ ' ' = \y{x) for all x G (-1, +1). 



The spectral properties of the equation ()4.2|1 are considered in the Hilbert function space 
L2(-l,+l). 

The van Hoeij method of searching for factors and solutions of this differential equation 
depends on the use of the computer algebra program Maple; see [121 ^^^1 [Tlj . 

Following the solutions of the fourth-order Bessel-type obtained in [HI Section 2] , and then 
in [Ij on using the methods of [H], it has proved possible to express all solutions of the 
equations ()4.1|) and ()4.2j) in terms of the confluent hypergeometric functions. 

The special functions involved are the Legendre functions Py{z) and Qy{z) as defined in 
the compendium Pfl Chapter 8, Section 8.1]. 

Define the functions P{v,x) and Q(z/, x), for the range of variables given in ^ Chapter 
8], by: 

(4.3) P(z/,x):=P°(x) 



(4.4) Q(z/,x):=Q°(z). 

Given the parameters A G (0, oo) and A G C define, noting the use of the ± symbol, 

(4.5) r^(A, A) := ±V4A2 - 4A + 1 + A and ^]^(A, A) := ^/5-8A + 4r±(A, A) 

We designate the four linearly independent solutions of ()4.1|1 . equivalently ()4.2j) . again the 
parameters A G (0, oo) and A G C, by 

(4.6) Ler{X, A; x) for all x G (-1, +1) and r = 1, 2, 3, 4. 
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These solutions are defined as follows, where d/dx denotes differentiation with respect to 
the independent variable x, 
(4.7) 

Lei(A, A] x) := 

-l{l + Q+{X,A))xP v^9 -8A + 4r+(A, A) + 4^]+(A, A) - \, a;) + 
[-(A + 3 - 4A + 4^2) + ^]+(A, A) - 3r+(A, A) + fi+(A, A)r+(A, A) + (A + 4A + 4A2)x2] x 



(4.8) 



(A + 4A + 4^2)-! Ap(i ^9 _ + 4r+(A, A) + 4fi+(A, A) - i, x) 



Le2(A,y4;a;) := 

1(1 + fi-(A,v4))a;P (^1^9 - + 4r-(A, A) + 4fi-(A, A) - i, a;) + 



[-(A + 3 - 4A + 4^2) + ^]-(A, A) - 3r-(A, A) + fi-(A, A)r-(A, A) + (A + 4A + 4A2)x2] x 

d 



(A + 4A + 4^2)-!— P(1^9 - + 4r-(A, A) + 4fi-(A, A) ~\,x). 



(4.9) 



Le3(A, A;x) := 

-i(l + ^]+(A, A))xg (1^9 - 8A + 4r+(A, A) + 4^]+(A, A) - + 

[-(A + 3 - 4A + 4^2) + fi+(A, A) - 3r+(A, A) + fi+(A, v4)r+(A, A) + (A + 4A + 4^2)^2] x 

d 



(4.10) 



(A + 4A + 4^2)-!— g(i^9 - 8A + 4r+(A, A) + A) - \ 

Le^^X, A; x) : 



-|(1 + n-{X,A))xQ (^1 v^9 -8A + 4r-(A, A) + 4fi-(A, A) - + 
[-(A + 3 - 4A + 4^2) + ^]-(A, A) - 3r-(A, A) + Q'{X, A)r'{X, A) + {X + AA + 4A'^)x^] x 
(A + 4A + 4^2)-! Aq(1 ^9 _ + 4r-(A, A) + 4^]-(A, A) - i, x). 

(1) It has been verified, using Maple 9, that Lei, Le2, Le^, Le^ are solutions of the dif- 
ferential equations ()4.ip and ()4.2j) for the domain x G (— 1,+1), and the parameter 
values A G (0, oo) and A G C; this verification follows from the original results of van 
Hoeij . 

(2) The factor (A + AA + 4^2)^^ which appears in the derivative term of all four of the 
solutions Lcr, for r = 1,2,3,4, presents an apparent singularity in these solutions 
when A = —{4A + 4^2); however Maple 9 determines that the other factors 

[-(A + 3 - 4A + 4^2) + Q^{X, A) - 3r±(A, A) + Q^{X, A)r^{X, A) + {X + 4A + 4A^)x^] 

both have a zero at this value A = —{AA + 4^42). 

5. JACOBI-TYPE DIFFERENTIAL EQUATION 

This note concerns information about the van Hoeij solutions to the fourth-order Jacobi- 
type linear ordinary differential equation. 
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The Jacobi-type differential equation was discovered by H.L. Krall, see JUl and ^1], and 
subsequently studied by other authors, see [Oj, [12] and [Ij. 

The differential equation may be written in two forms; here the parameters A G (0, oo) 
and a G (—1, oo), and the spectral parameter A G C: 

(1) The Frobenius form: 



(5.1) 



(1 - x2)2y(4)(a;) - 2(1 - x^){{a + 4)x + a)y^^\x) 
+ (1 + x) ((4A2" + + 9a + 14)x + (-4A2" + - 3a - 10)) y"{x) 
+ ((4Aa2° + 8A2" + 2a2 + 6a + 4)a; + (4Aa2" + 2a2 + 6a + 4)2/'(x) 

= Xy{x) for all a; G (—1, +1). 

(2) The Lagrange symmetric form (formally self-adjoint form): 

r J[y]{x) := ((1 - x)-+2(i + x)Y{x))" 

(5.2) < -((1 -2;)°+i((4A2" + 2a + 2)a; + 4A2° + 2a + 6)?/'(x))' 

[ = A(l -x)"?/(a;) for all X G (-1,+1). 

The spectral properties of the equation ()5.2|) are considered in the weighted Hilbert func- 
tion space L^((-l, +1); (1 - x)°). 

The van Hoeij method of searching for factors and solutions of this differential equation 
depends on the use of the computer algebra program Maple; see [13] and 

Following the solutions of the fourth-order Bessel-type obtained in [HI Section 2] , and then 
in [1] on using the methods of ^3], it has proved possible to express all solutions of the 
equations ()5.1|) and ()5.2|) in terms of hypergeometric functions. 

The special functions involved are the JacobiP functions; these functions are defined in 
terms of the hypergeometric functions pFq] for certain values of the parameters involved in 
this definition, they reduce to the classical Jacobi polynomials P^'^\ The notation JacobiP 
is not to be found in ^1 Chapter 22] but is given in the program MapleQ; here we shorten 
the notation JacobiP and define the function JP as follows, in terms of the classical P and 
2F1 functions: 

JacobiP(i^, a, /9; z) = JP{i', a, /5; z) 

(5.3) { P(z/ + a + l) 



P(z/+l)P(a + l 



-2Fi{-u, z/ + a + /? + l;a + l;(l- z)/2). 



Here, the independent variables z/, a, f3, z are in the complex field C, but for the results in 
this note we make the restrictions 

u eC, ae (-1, 00), /3 = 0, z e M. 

When z/ = n G No then we have the connection, see [H Section 22.2] and [3 Lecture 2, 
(2.2)], 

p('''f^\x) = JP{n,a,(3]x) for all a,/3 G (-l,oo) and all x G (-1,+1). 

The Jacobi-type orthogonal polynomials were discovered by H.L. Krall, see JU] and |llj . 
and later studied by A.M. Krall, T.H. Koornwinder and L.L. Littlejohn. These polynomials 

are denoted by P^j^{x) for a G (—1, 00), A G (0, 00), n G No and x G (—1, +1); each P"^^ is 
real-valued on (—1, -|-1) and of degree n. 
These polynomials have the properties: 
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(i) For all n G No the polynomial P"^ is a solution of the differential equation, see ()5.2|) 

(5.4) JiPnA"") = ^"(«' - ^TPnA^"") ^ ^ (-1' +1) 

where the eigenvalue A„(a,y4) is determined by 

(5.5) Xn{a, A) =n{n + a + l){n^ + {a + l)n + + a) for all n G Nq. 

(a) The collection {-P°^(-) : n G Nq} is a complete, orthogonal set in the Lebesgue- 
Stieltjes function space ((— 1, +1); /ijac) , where the Borel measure /xjac is deter- 
mined by the monotonic non- decreasing function /ijac, for all a G (— l,oo) and all 
A G (0, oo), 

Ajac(a;) := —1/2 for all x G (— oo, —1) 

(5.6) := ^|2"«- (1 tor all 

:= — r for all x G (+1, oo). 

( 2{a + l) 

The orthogonality of the set {PnA '■ ^ ^ -^o} then appears as 



^^^,A(-i)^n,A(-i) + 1^ l^'p:^jx)p:jx){i-xr dx 



'[-h+l] 

(5.7) 

2" 

for all m, n G No with m ^ n. 



The van Hoeij solutions for the general differential equation ()5.2j) begin with solutions to 
the set of special equations ()5.4|) when the general spectral parameter A takes the value of 
one of the eigenvalues {A„(a, A) : n E No}. 

Given n G No define, for all x G (—1, +1), with JP given by ()5.3p : 

(5.8) Si^n{x) := (na + 2A2" + n + n2)JP(n,a,0;x) 

+ (1 — x) — JP{n, a, 0; x). 
dx 

Then the van Hoeij methods determine that Si^n is a solution of ()5.1|) for this value of n, i.e. 
for A = \n{(y, A). 

Remark 5.1. From the definition ()5.8j) it would appear the this solution Si^n is linearly 
dependent upon the Jacob i-type polynomial P^a- 

To proceed further, given A G C, let {pr(A) G C : r = 1, 2, 3, 4} denote the roots of the 
quartic polynomial 

(5.9) p(p + a + l)(p2 + {a + l)p + 4A2° + a) - A = 0. 
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Maple gives explicit representations for these roots as follows: 

pi(a, A, A) 
P2{a,A,X) 
P3{a,A,X) 



(5.10) 



-a 



1 + v^a^ + I _ 8^2" + 2^{a, A, A) 
a - 1 - y/a"^ + 1 - 8A2° + 2^{a, A, A) 
a - H 



«2 + 1 _ - 2^(a, A A) 
P4(a, ^, A) = I (-a - 1 - _^ ^ _ 3^2" - 2^{a, A, A)) 



where 



^(a, A A) := + 8aA2" + 16^222" + 4A. 
The van Hoeij solutions to the Jacobi-type equation, see ()5.10|) . 
(5.11) J[y]i^) = A(l - x)"y(x) for all x G (-1, +1) 

now develop as follows. 

With A G C, replace n in the solution Si^n with the first root pi(A) taken from ()5.1()|1 : this 
yields a solution J7i in terms of the function JP, for all a G (—1, 00) and all A G (0, 00), 



(5.12) 



Ji{a, A, A; x) 



[a — ^(a, A, A)] JP{pi{a, A, A), a, 0; x) 

— (1 — x) — JP(pi(a, A, A), a, 0; x) for all x G (—1, +1). 
ax 



Similarly, for A G C, there is a solution J'2 derived from the third root P3(A), for all 
a G (—1, 00) and all A G (0, 00), 



(5.13) 



J2(a, A, A; x 



[a + ^(a. A, A)] JP^p^^a, A, A), a, 0; x) 

-(1 - x)-^JP{p3{a, A, A), a, 0; x) for all x G (-1, +1). 
ax 



Given n G Nq define, for all x G (—1, +1), there is a second solution 5'2,„, similar to Si^n 
of dHH), defined by 

52,n(x) := [((n + l)a + A2"+^ + n + n2)](x-l)-"JP(-l-n,-a,0;x) 

(5.14) ' / M / 

^ ' -{x -lY~'^—JP{-l-n,-a,0]x). 

ax 

Then the van Hoeij methods determine that S2,n is a solution of ()5.H) for this value of n, i.e. 
for A = An(a, A). 

With A G C, replace n in the solution S2,n with the first root pi(A) taken from ()5.10|1 : this 
yields a solution J'2, in terms of the function JP, for all a G (—1, 00) and all A G (0, 00), 
(5.15) 

J3{a,A,X;x) := --[a + ^(a, A, A)](-l + x)-"JP(-l - pi(a, A, A), -a, 0; x) 

+ (-1 + x)i-" — JP(-1 - pi(a, A, A), -a, 0; x) for all x G (-1, +1). 

fx J/ 
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Similarly, for A G C, there is a solution J74 derived from the third root P3(A), for all 
a G (—1, go) and all A G (0, oo), 
(5.16) 

j;(a,v4, A;x) := -^[a - ^(a, A, A)](-l + x)"" JP(-1 - ^3(0, A, A), -a, 0; x) 

d 

+ (-1 + a;)!-"— JP(-1 - p3(a, A, A), -a, 0; x) for all x G (-1, +1). 
dx 

6. Maple files 

The formal validity of these solution results were originally obtained by van Hoeij, following 
the methods given in and [14j, using the computer program Maple. These results have 
been subsequently confirmed by David Smith in the Maple program at the University of 
Birmingham. 

The Maple .mws files for the following three cases are named as follows: 

Jacobi-type-solutions.mws 
Laguerre-type-solutions.mws 
Legendre-type-solutions.mws 

Copies of these files may be obtained on application to David Smith at 

smithd@for.mat.bham.ac.uk 
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